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HIGHLIGHTS 


•  Preparation  and  characterization  of  nanostructured  rhodium  electrodes. 

•  Kinetic  study  of  the  hydrogen  oxidation  reaction. 

•  Correlation  of  experimental  results  through  a  unique  set  of  kinetic  parameters. 

•  Reaction  takes  place  under  both  Tafel-Volmer  and  Heyrovsky-Volmer  routes. 
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The  hydrogen  oxidation  reaction  was  studied  on  a  nanostructured  rhodium  electrode  at  different  rota¬ 
tion  rates.  The  electrode  was  prepared  via  sputtering  on  a  glassy  carbon  disc  support  and  it  was  char¬ 
acterized  by  X-ray  Photoelectron  Spectroscopy  (XPS),  Atomic  Force  Microscopy  (AFM)  and  cyclic 
voltammetry,  which  allowed  verifying  the  nanostructured  morphology  and  the  absence  of  any  phase 
other  than  metallic  rhodium.  The  real  surface  area  was  evaluated  by  CO  stripping  voltammetry.  Exper¬ 
imental  current  density  (j)  -  overpotential  (q)  curves  of  the  hydrogen  oxidation  reaction  were  obtained 
in  the  range  -0.015  V  <  q  <  0.25  V  at  different  rotation  rates  in  sulphuric  acid  solution.  They  were 
correlated  by  kinetic  expressions  derived  from  the  Tafel— Heyrovsky-Volmer  mechanism  and  thus  the 
kinetic  parameters  were  evaluated.  It  was  verified  that  over  this  overpotential  region  the  reaction  pro¬ 
ceeds  through  the  simultaneous  occurrence  of  the  Tafel— Volmer  and  the  Heyrovsky-Volmer  route. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  hydrogen  oxidation  reaction  ( hor )  has  received  considerable 
attention  in  relation  to  the  development  of  the  H2-O2  proton  ex¬ 
change  membrane  fuel  cell  (PEMFC)  [1].  Therefore,  the  elucidation 
of  the  factors  that  define  the  electrocatalytic  behaviour  for  this  re¬ 
action  is  essential  for  the  development  of  new  and  cheaper  elec¬ 
trode  materials.  Kinetic  studies  of  the  hydrogen  electrode  reaction 
have  been  carried  out  on  different  metals  preferably  in  the  cathodic 
overpotentials  range,  which  corresponds  to  the  hydrogen  evolution 
reaction  (her)  [2  ].  Starting  from  these  results,  it  has  been  established 
that  platinum  has  the  highest  electrocatalytic  activity,  as  it  can  be 
derived  from  its  position  on  the  top  of  the  volcano  curve,  which 
relates  the  exchange  current  density  (/„)  with  the  adsorption  free 
energy  of  the  reaction  intermediate  [3-6],  On  the  other  hand,  the 


*  Corresponding  author.  Tel:  +54  342  457  1164x2519;  fax:  +54  342  4536861. 
E-mail  address:  achialvo@fiq.unl.edu.ar  (AC.  Chialvo). 

0378-7753 /$  -  see  front  matter  ©  2014  Elsevier  B.V.  All  rights  reserved. 
http://dx.doi.org/10.1016/jjpowsour.2013.12.095 


number  of  kinetic  studies  of  the  hor  is  significantly  less,  most  of 
them  carried  out  on  Pt,  on  smooth  electrodes  [7—10]  as  well  as  on 
supported  nanoparticles  [11— 13],  Among  the  scarce  kinetic  data  of 
the  hor  on  other  noble  metals,  it  can  be  cited  the  evaluation  of  the 
kinetic  parameters  on  Pd  [14,15],  Ru  [16]  and  Ir  [17], 

There  are  only  a  few  studies  of  the  hydrogen  electrode  reaction 
on  rhodium.  The  first  antecedent  corresponds  to  an  evaluation  of 
the  equilibrium  polarization  resistance  [18],  Another  work 
measured  the  exchange  current  density  and  the  Tafel  slope  (b)  for 
the  hor  as  a  function  of  hydrogen  partial  pressure  [19],  These 
experimental  kinetic  parameters  were  also  measured  for  the  her  on 
acid  solution  [20].  More  recently  this  reaction  was  studied  by 
faradaic  impedance,  including  the  hydrogen  absorption  process 
[21],  Therefore,  there  has  been  no  report  about  kinetic  studies  of 
the  hydrogen  oxidation  reaction  on  Rh  electrodes.  On  this  context, 
the  present  work  deals  with  the  evaluation  of  the  elementary  ki¬ 
netic  parameters  of  the  Volmer-Heyrovsky-Tafel  mechanism  on 
nanostructured  Rh  electrodes  through  the  correlation  of  experi¬ 
mental  current  density  (j)  vs.  overpotential  (77)  curves,  measured  in 
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sulphuric  acid  solution  under  steady  state  and  controlled  mass 
transport  conditions. 

2.  Experimental  details 

2.1.  Electrode  preparation 

The  working  electrodes  were  prepared  via  sputtering  on  a  glassy 
carbon  (GC)  substrate  from  a  rhodium  target  provided  by  Good- 
fellow  Corp.  in  an  argon  atmosphere  (0.1  mbar),  using  a  sputter 
coater  Emitech  K500X,  operated  at  30  mA  during  4  min.  Previously 
to  Rh  deposition,  the  glassy  carbon  (SPI-Glas™  25)  substrate  was 
mirror  polished  and  subjected  to  ultrasonic  cleaning  in  ultrapure 
water  for  5  min  and  then  was  voltammetrically  characterized  to 
ensure  a  clean  and  reproducible  deposition  surface. 

2.2.  Electrode  characterization 

The  surface  morphology  of  the  Rh  electrodes  was  characterized 
by  Atomic  Force  Microscopy  (AFM).  AFM  images  were  obtained 
using  a  multi-technique  Agilent  microscope  model  5400,  operated 
in  tapping  mode  and  processed  with  the  software  WSxM  6.2. 

The  Rh  electrode  was  also  physicochemically  characterized  by 
X-ray  Photoelectron  Spectroscopy  (XPS)  using  a  multitechnique 
system  (SPECS)  equipped  with  a  dual  Mg/Al  X-ray  source  and  a 
hemispherical  PHOIBOS  150  analyzer  operating  in  the  fixed 
analyzer  transmission  (FAT)  mode.  The  spectra  were  obtained  with 
pass  energy  of  30  eV;  a  Mg  Ka  X-ray  source  was  operated  at  200  W 
and  12  kV.  The  working  pressure  in  the  analysing  chamber  was  less 
than  5.9  x  10-7  Pa.  The  XPS  analyses  were  performed  in  the 
spectral  region  corresponding  to  Rh  3d  core  level.  The  data  treat¬ 
ment  was  performed  with  the  Casa  XPS  program  (Casa  Software 
Ltd.,  UK). 

2.3.  Electrochemical  measurements 

All  the  electrochemical  measurements  were  performed  in  a 
three  electrodes  cell  specially  built  for  the  use  of  a  rotating  elec¬ 
trode  and  with  a  particular  design  of  the  gas  saturator  at  298  K.  The 
electrolytic  solution  was  0.5  M  H2SO4,  prepared  with  ultra-pure 
water  (PureLab,  Elga  LabWater).  The  working  electrode  was 
mounted  in  a  teflon  holder,  with  an  exposed  geometric  area  of 
0.07  cm2.  It  was  rotated  through  the  use  of  a  rotating  disk  Radi¬ 
ometer  EDI  10K  system.  Furthermore,  a  large  area  platinum  wire 
acted  as  counter  electrode  and  a  hydrogen  electrode  in  the  same 
solution  (RHE)  was  used  as  reference  electrode.  The  working  elec¬ 
trode  was  electrochemically  characterized  by  cyclic  voltammetry 
using  a  potentiostat-galvanostat  Wenking  POS2  controlled  by  an 
interface  Advantech  PI1710HG  and  the  software  Labview  8.  The 
voltammograms  were  obtained  at  0.1  V  s-1  in  nitrogen  saturated 
solution  between  0.0  and  1.3  V. 

The  electrocatalytic  activity  towards  the  hydrogen  oxidation 
reaction  was  evaluated  through  the  current  (!)  —  overpotential  (rj) 
response  in  hydrogen  saturated  solution.  The  polarization  curves 
were  obtained  at  different  rotation  rates  in  the  range  900  <  w/rpm 
<4900.  A  potential  program  was  applied,  which  initiated  with  a  3  s 
step  at  0.0  V,  followed  by  a  5  s  step  to  the  desired  overpotential 
value.  In  this  last  period,  readings  of  the  current  value  were  made 
each  0.1  s  and  the  mean  value  of  the  current  data  measured  in  the 
last  2  s  was  assigned  to  the  step  overpotential.  Then  the  program 
was  repeated  for  each  p  value,  which  was  varied  in  the 
range  -0.015  <  77/V  <  0.25.  The  electrolyte  solution  was  renewed 
after  each  experiment. 

Finally,  the  real  area  of  the  electrode  was  determined  by  CO 
stripping  voltammetry.  CO  adsorption  was  carried  out  in  an 


auxiliary  cell  in  a  CO  saturated  solution,  holding  the  electrode  po¬ 
tential  at  0.05  V  during  15  min.  Then,  the  voltammetric  stripping 
was  carried  out  at  0.05  V  s_1  in  nitrogen  saturated  solution  be¬ 
tween  0.0  and  1.3  V. 

3.  Results  and  discussion 

3.1.  Morphological  and  physicochemical  characterization 

The  surface  morphology  was  analysed  by  AFM  (Fig.  1 ).  It  can  be 
observed  that  the  Rh  layer,  deposited  by  sputtering  on  a  glassy 
carbon  substrate,  consists  of  uniformly  distributed  nanoparticles.  In 
order  to  estimate  the  particle  size,  the  height  profile  was  analysed 
in  different  places  on  the  electrode  surface.  Two  of  them  are 
illustrated  in  the  inset  of  Fig.  1.  The  average  height  of  the  nano¬ 
particles  obtained  from  these  profiles  is  ranging  between  2  and 
6  nm. 

The  physicochemical  characterization  of  the  Rh  electrode  was 
carried  out  by  XPS  analysis.  Fig.  2  illustrates  the  spectrum  in  the 
range  of  binding  energy  comprised  between  303  and  318  eV 
(continuous  line).  It  can  be  observed  two  peaks  located  at  307.0  eV 
and  311.8  eV  respectively,  which  are  in  agreement  with  the  values 
reported  for  3ds/2  and  3d3/2  signals  of  metal  rhodium  [22],  The 
deconvolution  of  these  peaks  is  also  shown  (dotted  lines). 

3.2.  Electrochemical  characterization  and  real  area  evaluation 

Fig.  3  shows  the  voltammetric  profile  of  the  nanostructured  Rh 
electrode  obtained  in  0.5  M  H2SO4  solution  at  0.05  V  s-1  under 
nitrogen  bubbling.  The  electrode  was  previously  subjected  to  a  CO 
adsorption  carried  out  in  an  auxiliary  cell  in  a  CO  saturated  solution 
at  0.05  V  during  15  min.  The  first  cycle  corresponds  to  the  stripping 
of  the  adsorbed  CO,  where  it  can  be  observed  a  sharp  peak  at 
approximately  0.725  V.  The  complete  electrooxidation  of  the 
adsorbed  CO  is  produced  in  this  first  anodic  sweep,  as  it  can  be 
verified  through  the  second  voltammetric  cycle,  which  reproduces 
the  blank  voltammogram  of  metallic  rhodium  in  this  electrolyte 


Fig.  1.  AFM  image  of  the  nanostructured  Rh  electrode.  Insets:  height  profiles. 
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Binding  Energy  /  eV 


Fig.  2.  XPS  Rh  3d  spectra  of  the  nanostructured  Rh  electrode  (continuous  line)  and 
deconvolution  of  the  peaks  (dotted  lines). 


solution.  It  can  be  appreciated  that  the  voltammetric  profile,  and  in 
particular  the  potential  region  where  the  hydrogen  adsorption 
takes  place,  is  very  similar  to  those  found  in  the  literature  [23,24], 
Moreover,  this  result  is  in  agreement  with  the  XPS  analysis,  where 
only  metallic  rhodium  was  found. 

The  voltammetric  charge  corresponding  to  the  electro¬ 
oxidation  of  the  adsorbed  CO  was  used  for  the  determination  of 
the  real  electrode  area.  This  charge  was  evaluated  as  the  differ¬ 
ence  between  the  first  and  second  cycle  in  the  potential  range 
0.40  <  E/V  <  1.3,  resulting  a  value  equal  to  110.5  pC.  Taking  into 
account  that  for  Pt  [25],  Pd  [26]  and  Ir  [17],  the  surface  coverage 
corresponding  to  the  CO  saturation  is  66%  of  the  total  surface 
sites,  it  should  be  reasonable  to  consider  that  this  value  is  also 
applicable  for  the  case  of  Rh.  On  this  basis,  the  voltammetric 
charge  of  one  CO  monolayer  is  equal  to  287.7  pC  cm-2  and 


in  the  range  -0.015  V  <  ij  <  0.25  V  at  different  rotation  rates  (indicated  in  the  figure). 
0.5  M  H2S04;  298  K. 


electrode.  0.5  M  H2S04;  298  K.  to  =  (A)  900;  (•)  1600;  (■)  2500;  (★)  3600;  (♦) 
4900  rpm. 


Fig.  3.  Voltammogram  of  the  Rh  electrode  ai 
Second  cycle.  0.05  V  s”1;  0.5  M  H2S04,  298  K. 


CO  adsorption,  (a)  First  cycle;  (b) 
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therefore  the  real  area  of  the  Rh  electrode  is  0.384  cm2.  More¬ 
over,  as  the  reaction  under  study  has  a  strong  diffusion  contri¬ 
bution,  which  depends  on  the  geometric  area,  meanwhile  the 
kinetic  reaction  takes  place  on  the  real  surface  area,  the 
relationship  between  both  values  is  an  important  parameter  for 
the  correct  analysis  of  the  reaction.  This  relationship,  called 
active  area  factor  (/oa),  is  in  this  case  equal  to  5.486.  It  appears  in 
the  kinetic  equations  for  the  hydrogen  oxidation  reaction  (item 
3.3). 

3.3.  Hydrogen  oxidation  reaction 

The  evaluation  of  the  experimental  dependence  current  (/)  vs. 
7j  for  the  hydrogen  oxidation  reaction  on  a  nanostructured  Rh 
electrode  was  carried  out  in  0.5  M  H2SO4  solution  under  hydrogen 
bubbling  at  different  rotation  rates  (900  <  w/rpm  <4900).  After 
the  open  circuit  potential  reached  the  equilibrium  value 
(0.0  ±  0.0004  V  vs.  RHE),  a  potential  program  was  applied,  which 
initiated  with  a  3  s  step  at  0.0  V,  followed  by  a  5  s  step  to  the 
desired  overpotential  value.  Then  the  program  was  repeated  for 
each  t]  value,  which  was  varied  in  the  range  -0.015  <  77/V  <  0.25. 
The  current  on  time  (t)  response  of  the  Rh  electrode  resulting 
from  the  application  of  the  potential  program  is  shown  in  Fig.  4, 
for  all  the  rotation  rates  analysed  (900,  1600,  2500,  3600  and 
4900  rpm).  It  can  be  observed  that  current  quickly  reached  the 
steady  state  condition.  Readings  of  the  current  value  were  made 
each  0.1  s  and  the  mean  value  of  the  current  data  measured  in  the 
last  2  s  was  assigned  to  the  step  overpotential.  Starting  from  these 
results  and  taking  into  account  the  real  electrode  area,  the  cor¬ 
responding  j(i))  dependences  on  steady  state  were  evaluated  for 
each  rotation  rate,  which  are  shown  in  Fig.  5  (symbols).  These 
experimental  curves  were  correlated  with  the  system  of  equa¬ 
tions  resulting  from  the  resolution  of  the  kinetic  mechanism  of 
Tafel-Heyrovsky-Volmer  on  steady  state,  which  is  valid  for  the 
hydrogen  electrode  reaction  on  a  metallic  surface.  The  corre¬ 
sponding  steps  can  be  written  as, 

H2+2Sf±2Had  Tafel  (la) 

H2  +S<=sHad  +H+  +  e~  Heyrovsky  (lb) 

Had  <=s  H+  +  e~  +  S  Volmer  (lc) 


The  implicit  equation  for  the  variation  of  the  surface  coverage 
(0)  of  the  reaction  intermediate  on  overpotential  is  given  by, 


peU(J/).-(l(~yj  +vi1e-u(»-n^ 
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In  equations  (2)  and  (3),  V;  is  the  reaction  rate  of  the  step  i  (i  =  T, 
H,  V),  a i  (i  =  V,  H )  is  the  symmetry  factor  of  the  step  i,  A  is  the 
symmetry  factor  of  adsorption,  u  (in  RT  units)  is  the  energy  of 
interaction  between  the  adsorbed  hydrogen  atoms  and  B  is  the 
Levich  constant.  Superscript  e  indicates  equilibrium  and  /  =  F/RT. 
The  experimental  j{t])  curves  for  each  w  value  can  be  correlated 
using  equations  (2)  and  (3),  applying  a  non  linear  least  squares 
regression  method.  However,  it  is  necessary  first  to  evaluate  the 
value  of  constant  B.  The  relationship /afl/B  can  be  evaluated  from  the 
values  of  the  maximum  current  (j max)  observed  for  t]  >  0.20  V,  as  it 
is  related  to  the  rotation  rate  as  follows  [28], 


jmax(w)  wl/2  (4) 

where  is  the  limiting  kinetic  current  density  of  the  Tafel  step, 
which  represents  the  mixed  behaviour  of  the  hydrogen  oxidation, 
with  simultaneous  kinetic  and  diffusion  contributions.  It  is  given  by 
the  following  expression  [28], 


kin  _2 Fvfe2^ 

•W"  (i_rf 


(5) 


where  S  represents  a  site  on  the  electrode  surface  in  which  the 
reaction  intermediate  Had  can  be  adsorbed.  The  experimental 
dependences  j{ri,  w)  were  described  with  the  kinetic  expression 
resulting  from  the  rigorous  resolution  of  the  kinetic  mechanism 
on  steady  state,  considering  a  Frumkin  type  adsorption  of  the 
reaction  intermediate.  Details  of  the  resolution,  which  considers 
the  simultaneous  occurrence  of  the  two  independent  routes, 
Tafel-Volmer  and  Heyrovsky-Volmer,  can  be  found  elsewhere 
[[27]  and  references  cited  there].  Three  equivalent  expressions  for 
j(r),  w)  are  derived,  referred  to  the  real  or  active  area,  resulting 
from  the  mass  balance  in  steady  state.  The  equation  used  in  the 
present  work  was, 


Therefore,  the  constants /aa/B  and  jj^Jx  can  be  obtained  from  the 
slope  and  the  origin  ordinate  of  the  experimental  plot  vs.  w-1^2. 

Fig.  6  illustrates  such  plot,  where  it  can  be  observed  that  the 
experimental  points  (dots)  satisfy  the  required  linearity.  From 
linear  regression  (line)  the  following  values  were  obtained:  faa/ 
B  =  99.129  mA1  cm2  rpm1/2  and  (JmSx)-1  =  0.02707  mA  1  cm2.  The 
following  relationship  between  the  parameters  vf.  and  8e  is  also 
obtained  from  equation  (5), 

vf.[mol  cm-vij 1.914  10  7(1  -  0e)V2u(f  (6) 

Then,  the  experimental  dependences  j(?j,  w)  were  correlated 
using  equations  (2),  (3)  and  (6),  applying  a  non  linear  least  squares 
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regression  method  carried  out  with  the  software  Micromath  Sci¬ 
entific  3.0,  considering  the  symmetry  factors  av  =  uh  =  A  =  0.5.  The 
kinetic  parameters  v^veH,de  and  u  were  evaluated  from  the  correla¬ 
tion  and  the  value  of  vf  was  obtained  from  equation  (6).  Each  po¬ 
larization  curve  corresponding  to  a  particular  rotation  rate  was 
correlated  independently.  The  minimization  of  the  mean  squared 
error  led  for  all  the  rotation  rates  to  the  cancellation  of  the  inter¬ 
action  parameter  u,  which  means  that  the  adsorption  of  the  reac¬ 
tion  intermediate  follows  a  Langmuir  isotherm.  The  values  of  the 
kinetic  parameters  Py.vf]  and  9e  as  well  as  the  corresponding  mean 
values  are  shown  in  Table  l.The  simulations  of  the  dependences  j(rj, 
w)  were  done  with  this  unique  set  of  parameters,  which  are  shown 
in  Fig.  5  (lines).  It  can  be  appreciated  that  there  is  a  good  agreement 
between  the  experimental  and  simulated  curves  and  that  in  the 
analysed  potential  range  the  hydrogen  oxidation  is  verified  through 
both,  Tafel-Volmer  and  Heyrovsky— Volmer  routes,  being  the  two 
contributions  significant. 


literature  (0.316  mA  cm-2)  [30],  It  is  also  almost  two  orders  of 
magnitude  lower  than  that  corresponding  to  Pt  nanoparticles 
supported  on  GC  (45.9  mA  cm-2)  [11], 

Moreover,  the  equilibrium  polarization  resistance  (R£)  can  be 
also  evaluated  from  the  kinetic  parameters  [29], 

RT  (  44  +  vf,  +  4 

p  4F2  ^44  +  44  +  44 

The  resulting  value  is  R°  =  38.6  Q  cm2,  which  is  higher  than  that 
reported  byJ.P.  Hoare  and  S.  Schuldiner  (25.0  Q  cm2)  [18],  However, 
it  can  be  appreciated  in  Fig.  1  of  such  reference  that  the  dependence 
j( 7))  is  not  linear  and  that  from  the  tangent  line  at  the  origin  the 
value  Rp  —  44.0  Q  cm2  is  obtained,  which  is  more  in  agreement  with 
the  present  result.  As  it  is  expected,  this  value  is  also  significantly 
higher  than  that  on  Pt  (0.355  Q  cm2)  [11], 

4.  Conclusions 

The  hydrogen  oxidation  reaction  was  studied  on  a  rotating  disc 
electrode  of  nanostructured  rhodium  supported  on  glassy  carbon, 
being  the  electroactive  area  evaluated  by  CO  stripping  voltamme¬ 
try.  The  kinetic  parameters  resulting  from  the  rigorous  resolution  of 
the  Tafel-Heyrovsky-Volmer  mechanism  were  evaluated  from  the 
correlation  of  the  experimental  current  density— overpotential 
curves  measured  at  different  rotation  rates.  It  was  demonstrated 
that  in  the  analysed  overpotential  region  the  reaction  takes  place 
through  the  Tafel-Volmer  route  and  the  Heyrovsky-Volmer  route. 
It  was  found  that  both  contributions  are  significant  and  that  there  is 
no  prevalence  of  one  route  over  the  other.  Self  consistent  values  of 
the  exchange  current  density  and  the  equilibrium  polarization 
resistance  of  the  hydrogen  electrode  reaction  were  also  evaluated 
from  these  kinetic  parameters.  These  results  were  discussed  and 
compared  with  values  reported  in  the  literature. 

Acknowledgements 

The  authors  wish  to  acknowledge  the  financial  support  received 
from  ANPCyT,  CONICET  and  UNL.  Thanks  are  also  given  to  ANPCyT 
for  the  purchase  of  the  SPECS  multitechnique  analysis  instrument 
(PME2003-8). 


3.4.  Electrocatalytic  activity 
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